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Boero ef al. have recently proposed jy a model for 
the E' center in quartz and silica, based on ab initio 
studies of the singly positive charge state of the oxygen 
vacancy Vo- They find that in the +1 charge state, Vo 
undergoes a distortion to a puckered configuration whose 
calculated hyperfine activity is consistent with the ob- 
served data. The purpose of this Comment is to point 
out that recent ab initio density-functional calculations 
on the oxygen vacancy in a— quartz [Q rule out the pos- 
sibility of identifying the +1 charge state of the vacancy 
with the E ' center as suggested in |TJ . Our results [Q in- 
dicate the Q = —3 charge state of Vo as another possible 
candidate. 
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FIG. 1. Calculated formation energies for the various 
charge states of the oxygen vacancy in a— quartz. 

The formation energies of the oxygen vacancy are dis- 
played in Figjl as a function of the Fermi level (whose 
zero is chosen at the valence band top). The oxygen va- 
cancy is a negative-£/ center whose ground state can have 
the charge states Q= [+3,0,-2,-3]. The singly charged 
states, in particular the +1 state, are never stable against 
capture or release of further electrons. Therefore the 
singly positive vacancy cannot be invoked as a candi- 
date E ' center. This remains true accounting for the 
puckering energy gain of 0.3 eV B. 

Our calculations are technically quite comparable to 
those in Qj. We studied a single oxygen vacancy in 36- 
atom a— quartz supercells using conjugate-gradient total 
energy minimization with ultrasoft pseudopotentials ||, 
a plane-waves basis cut off at 20 Ryd, and one special 
k-point. No symmetry restriction is imposed, and all 
geometries are fully relaxed. The theoretical lattice con- 
stants (deviating less than 1% from experiment) are used. 
The energy zeroes in the charged supercells are carefully 
aligned pj. 

The Q = — 3 state of Vo is characterized by a ma- 


jor lattice distortion involving three vacancy-neighboring 
tetrahedral units (Fig.||). While the Q = and Q = +3 
states fail to match most of the magnetic or optical prop- 
erties of the E ' center, the Q — —3 state satisfies most re- 
quirements for identification with E'. Calculated optical 
transition energies agree well with absorption and emis- 
sion bands commonly associated with E'. The main fea- 
ture of the charge density of the unpaired electron in the 
defect state (not shown due to space constraints) is a pro- 
nounced lobe associated with atom S13 (see FigJ|), and 
pointing approximately towards Sii . The angle formed 
by this dangling lobe with the neighboring O atoms is 
^100° in accordance with experimental estimates. Two 
additional, weaker density features are associated with 
Sii and Si2- The overall picture seems compatible with 
the observation of one strong, and two weak hyperfine sig- 
nals. Calculated hyperfine parameters are unfortunately 
not available at the moment. 



FIG. 2. Equilibrium configuration of Vo for Q — — 3 
(thick lines, black atoms) and Q = (thin lines, grey atoms). 

In summary the +1 state of the oxygen vacancy in 
quartz cannot be associated with the E' center. An al- 
ternative candidate is the -3 charge state of the same 
center. 
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